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We investigate low-lying octupole response of neutron-rich Ni isotopes beyond the N=50 
J~ . shell closure using the Skyrme-Hartree-Fock-Bogoliubov mean-fields and the continuum 

• quasi-particle random phase approximation. Performing detailed numerical analyses em- 

QQ ' ploying the Skyrme parameter set SLy4 and a density-dependent delta interaction of the 

mixed type, we show that a neutron mode emerges above the neutron separation energy 
as a consequence of the weak binding of neutrons and it exhibits strong influences of the 
di-neutron correlation. 

^ ! §1. Introduction 

The di-neutron correlation, a spatial pair correlation with small correlation 
length among neutrons, has been one of the central themes of the physics of two- 
neutron halo nuclei such as ^^Li.^^"^^ Although an affirmative experimental signature 
of the neutron spatial correlation in ^^Li is obtained rather recently,®-* many theo- 
, retical predictions have been accumulated, concerning not only two-neutron halo 

nuclei^^"^-*'^^"^^-* but also other neutron-rich systems including the surface area of 
■ medium- and heavy- mass neutron rich nuclei^^^' and dilute neutron matter. ^^^"^^^ 

CD I Indeed the analysis of the spatial structure of the neutron Cooper pair in dilute 

matter^'')"^^) has revealed a mechanism that the di-neutron correlation may emerge 
generically. Namely the spatial di-neutron correlation originates from a strong cou- 
pling feature of the neutron pair correlation, which becomes significant at low den- 
^ ' sities because of the strong momentum dependence of the attractive nuclear force in 

channel. Furthermore the induced pairing interaction caused by the exchange of 
surface phonons is claimed to give additional contribution to possible enhancement 
of the spatial correlation in finite nuclei.^'''' 

The influence of the di-neutron correlation on nuclear structure could be mul- 
tifold. Namely it may influence not only the ground state but also various modes 
of excitations and dynamics. Attentions, however, have been focused so far mostly 
on soft dipole excitation in two-neutron halo nuclei, ^^"^^'^^^ and rather recently, in 
medium- and heavy-mass neutron rich nuclei. ^^^'^^^ Concerning the latter case, we 
have shown in our previous study^^^ that strong influence of the di-neutron corre- 
lation in the soft dipole excitation of proton semi-magic neutron-rich Ca and Ni 
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nuclei. Considering the possible generality of the di-neutron correlation, it is ex- 
pected that this correlation may also emerge in other multipole modes of excitation. 
We anticipate it in multipole excitations with negative and natural parity because 
the di-neutron correlation involves a configuration mixing among single-particle or- 
bits with different parities. In the present paper, therefore, we examine low-lying 
octupole modes of excitation. 

There exists only a few investigations of low- lying octupole excitations in neutron- 
rich nuclei while in stable nuclei the low-lying 3^ state with a character of octupole 
surface vibration is well established. ^^^'^^^ Experimental properties of the low-lying 
octupole modes in neutron-rich nuclei are known very little, with exceptions such as 
20q 23) ^]^g theoretical side, it is predicted^^^"^^^ that neutron halo nuclei exhibit 
a neutron mode with large octupole strength in the continuum region above the neu- 
tron threshold, resulting from transitions from weakly bound orbits to continuum 
orbits, and thus having a completely different character with the surface vibration. 
Analyses based on the random phase approximation (RPA),^^)"^*^) which can de- 
scribe both the continuum states and the collectivity, have shown coexistence of the 
collective surface vibration and neutron continuum strength near the threshold in 
medium mass neutron rich nuclei. Note however that models of Refs.^^^"^*^^ do not 
include the pair correlations, and analyses are limited to nuclei with doubly closed 
shell configurations, such as ^'^Ca, ^^O and ^^'^^Ni. The quasi-particle random phase 
approximation including the pair correlation effect is employed in Ref.^^^ to describe 
oxygen isotopes ^*^^^^0. There, however, the Hartree-Fock-|-BCS approximation is 
adopted instead of the Hartrcc-Fock-Bogoliubov scheme, and effects of the pair cor- 
relation are not discussed intensively. In contrast to these preceding works, we would 
like to focus in the present paper on roles of the pair correlation and the possibility 
of the di-neutron correlation especially in the low-lying octupole correlation. For 
this purpose we use the continuum quasi-particle random phase approximation (the 
continuum 

QRPA)i5).i8)>3i) which can describe both the continuum, the pairing, 
and the collectivities in neutron rich medium-mass nuclei. We choose nickel isotopes 
80-86]\jj for numerical analysis, and perform detailed investigation of ^^Ni which is 
chosen as an representative example. A preliminary report of this work is seen in 
Ref.32) 

§2. Skyrme-HFB plus continuum QRPA method 

In the present study we assume the spherical symmetry of the ground state as 
we analyze semi-magic nickel isotopes. We first construct the spherically symmetric 
ground state and associated self-consistent mean-fields by means of the coordinate- 
space Skyrme-Hartree-Fock-Bogoliubov (Skyrme-HFB) method, ^^^'^^^ with which we 
can describe properly the spatially extended wave functions of weakly bound quasi- 
particle orbits. We employ the Skyrme-force parameter set SLy4,^^) which has been 
extensively used for studies of the neutron-rich nuclei. We shall also use an- 
other parameter set SkM**^"^)'^^) for comparison. In our previous study of the dipole 
excitation^^) a Woods-Saxon potential was adopted, but here we perform fully self- 
consistent HFB calculations. Concerning the effective pair interaction, we use the 
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density-dependent delta-type interaction (DDDI)^)'^^-*'^^) which is given by 

vpairir, r') = 1fo(1 - Pa) (l " V^) - r'). (2-1) 

We choose the so-called mixed type pairing, i.e. ry = 1/2, = 0.16 fm^^, concern- 
ing the density dependence, and we determine the strength of the pair interaction 
Vq in the same way adopted in Ref.^^^ Namely we determined Vq = — 315MeV fm'^ 
for nickel isotopes so that the average pairing gap for neutrons is in overall agree- 
ment with the experimental odd-even mass difference evaluated with the three-point 
formula'^s) (cf. FigH]). 

We describe excitation modes by means of the continuum QRPA. The contin- 
uum QRPA formalism is essentially the same as those in Refs.^^^'^^^ except that 
we here formulate it on the basis of the Skyrme-Hartree-Fock-Bogoliubov mean- 
fields and different residual interactions. Concerning the residual interaction in 
the pairing channel, we adopt the one derived from the functional derivatives of 
the pairing energy functional, as was in our previous study. ^^-^ Concerning the 
residual interaction in the particle-hole channel, we use the same Skyrme inter- 
action employed in constructing the ground state, but we take a Landau-Migdal 
approximation to the Skyrme interaction^^^"^^^ in order to make the systematic 
numerical calculations feasible. In practice, the particle-hole interaction is given 
by Vphir - r') = {{Fo/No){r) + {F^/No){r)T ■ T')}d{r - r') where Fq, F^ and iVo 
are the Landau-Migdal parameters'^^^"^^-' and the associated normalization factor, 
evaluated for the Skyrme interaction. Their expressions are given, for example, in 
Refs.^^^'^^^ The Landau-Migdal parameters are defined usually for symmetric nu- 
clear matter, but we treat the Fermi momentum kp in the expressions as a local 
quantity kp{r) = (37r^p(r)/2)^''^ related to the nucleon density p{r).^^^ Thus the 
force strength is density dependent, and hence depends on the spatial coordinate r. 

Given the HFB mean- fields and the residual interactions, we solve the linear 
response equation 

( 5p,L{r,u:) \ ( \ f j:,'^liir')6p,,Ur',u)/r'^+vll\r') 

5p+,gLir,u;) \= dr'\ Rf^^ir,r' ,u;) Kpa^r{r')6p+,,L{r' ,io)/r'^ 

V 6p.,,L{r,u;) y V / V -^pair{r')Sp^,,L{r' ,u;)/r'^ 

(2-2) 

to describe the correlated multipole response of a nucleus against an external field. 
Here Spqiir, uj), 6p+^qL{r, uj), and 6p-^qL{r, u>) are responses in the normal and abnor- 
mal densities with multi-polarity L, and Uph and Kpair are related to the interaction 
strengths (-Fo/-^o)(^')) iPo/^o)i^) ^(1 ~ VPi^)/ Po)- Note that we construct the 
response function R^^qpir, r', u>) in terms of products of two single-quasi-particle HFB 
Green's function summed over not only the discrete quasi-particle states but also the 
continuum quasi-particle states. This is possible because we use the exact HFB 
Green's function^^^ consisting of the regular and out-going wave solutions of the HFB 
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equation. We consider the external field 

y-*(r.) = ^rfy3o(fO and Y.l^rfY,o{r,), (2-3) 

i i 

for the isoscalar, neutron and proton strength functions dB[\L) / dE = B{XL; Ogs — > 
2>~)5{E—Ei) [XL =IS3, nS and p3, respectively), which can be evaluated as dB{XL)/dE = 
-^Im I T^gdrVext,qir)SpqL{r,uj) using the solution 5pqL{r,uj) of Eq.(l22]). 

Numerical details are as follows. The HFB equation is solved using the radial 
coordinate in a spherical box r = — r^ax- The Skyrme-HFB code is our orig- 
inal one, whose numerical procedures essentially follow those of Ref.^^-* We solve 
the radial HFB equation using the Runge-Kutta method instead of the Numerov 
method adopted in Refs.'^^^'^^-' since we need a consistent evaluation of derivatives 
of the quasi-particle wave functions in constructing the Green's function used in the 
continuum QRPA calculation. We have checked that results of our Skyrme-HFB 
code agree with those produced by the code HFBRAD.^^^ The continuum QRPA 
part is based on our previous version^^) employing the Woods-Saxon potential, but 
here we replace the Woods-Saxon potential by the Skyrme-HFB mean-fields. Also 
we implement the Landau-Migdal approximation of the Skyrme interaction as the 
residual interaction in the particle-hole channel. As the box size we choose r^ax = 
22 fm, and an equidistant discretization with Ar = 0.2 fm. Concerning the cut-off 
of the quasi-particle orbits, we set Imax = 17^ for the single-particle partial waves Ij, 
and Emax = 60 MeV with respect to the quasi-particle energies. The cut-off energy 
is a standard choice adopted in many Skyrme HFB calculations,^^^'^^^'^^^'^'^^''^^^ but 
slightly larger than that adopted in our previous Woods-Saxon calculation.^^) Con- 
cerning the orbital angular momentum cut-off, we find that the ground state and 
the octupole strength function have good convergence already at I = 12h, but to 
evaluate the transition densities we need larger Vs. Because we adopt the Landau- 
Migdal approximation, the consistency between the ground state and the excited 
states is partly broken. In order to minimize effects of the self-consistency breaking, 
we renormalize the strength of the particle-hole residual interaction as Vph — > / x Vph 
by a numerical factor / so that the spurious center-of-mass mode in the isoscalar 
dipole excitation has zero excitation energy. 

§3. Numerical analysis 

3.1. Octupole strength functions in ^"^Ni 

We shall first discuss a representative example ^^Ni in order to clarify basic 
features of low-lying octupole excitation. 

In Fig. [21 we show the isoscalar octupole strength in ^^Ni. Here the strength 
function is calculated with smoothing parameter e = 0.2 MeV, which is introduced 
in the linear response equation as the imaginary part of the excitation frequency 
u + ie. This means that the strength function is folded with an Lorentzian function 
with FWHM = 2e = 0.4 MeV. It is seen that there are essentially two groups of 
strength distribution; one around = 26 - 32 MeV and the low-lying distributions 
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below ~ 10 MeV. We regard the high energy group as the 3-huj high-frequency 
vibrational mode.^^^ In the low-energy group, the sharp peak at S ~ 4.15 MeV 
is most prominent, but another structure is also seen. We focus on the low-energy 
group in the following discussion. 

In Fig. [3] we show magnification of the isoscalar strength function in a low- lying 
region E = — 8 MeV, together with the neutron and proton octupole strength func- 
tions. The smoothing parameter is chosen to a smaller value e = 0.05 MeV in this 
case (and also in most of the following calculations unless mentioned explicitly). It is 
now clear in Fig. [3] that the low- lying strength consists of two structures. Apart from 
the sharp peak at E = 4.16 MeV, there exists a broad bump which emerges above 
the one-neutron separation energy (the one- and two-neuron separation energies are 
Ein,2n = 1-86 and 2.37 MeV as indicated with arrows). Although the two structures 
overlap in the same energy region, their characters are clearly distinguished as fol- 
lows. Firstly, the broad bump does not form a well defined peak, and we consider 
it as a kind of continuum mode where fast neutron decay takes place. It carries 
essentially only neutron strength. Secondly, the sharp peak has a small width even 
though it is located above the neutron separation energy. It can be regarded as a 
narrow resonance. It is also distinguished from the neutron mode by the fact that it 
carries sizable proton strength. 

In order to give more characterizations, we examined how the residual interac- 
tions influence these modes. Namely we performed three calculations where either or 
both of the particle-hole and pairing residual interactions are switched off. Results 
are shown in Fig. SI The sharp peak disappears when the particle-hole interaction 
is neglected, and it is influenced rather weakly by the pairing residual interaction. 
Thus the main origin of the sharp-peak mode is correlation due to the particle-hole 
residual interaction. It suggests that this mode could be the surface vibrational mode 
consisting of low-energy 1-tkj particle-hole transitions.^^-* On the contrary, the broad 
neutron mode exists even when the particle-hole residual interaction is neglected. It 
has different origin. 

It is useful to evaluate integrated sums of octupole strengths associated with the 
two modes. It is, however, not easy to evaluate them separately since the sharp-peak 
mode and the broad neutron mode overlap in the same energy region. Nevertheless 
we estimate them in the following manner. Concerning the broad neutron mode 
we integrate the strength functions in an energy region above the one-neutron sep- 
aration energy Ein < E < Ein + 1.5 MeV with an interval of 1.5 MeV, where the 
neutron strength dominates and the collective vibrational mode barely overlaps. The 
boundary energies are Ein = 1-86 and Ein + 1-5 = 3.36 MeV in the case of ®^Ni. 
Since the choice of the energy interval 1.5 MeV is rather arbitrary and probably 
small to cover the whole strength of the broad neutron mode, we expect that there 
may be an underestimate by a factor of up to about 2. Concerning the strength of 
the surface vibrational mode, we define an energy interval Ei < E < E2 where the 
strengths are integrated by noticing that this mode carries a proton strength, which 
is a character clearly distinguishable from the broad neutron mode. In practice, we 
define the boundaries Ei and E2 at which the proton strength function is 1.0% of 
the value at the peak energy Ep^ak- In ^^Ni, Ei and E2 are 3.50 and 5.00 MeV, 
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respectively, and there is no overlap between the two energy intervals. In Table [U we 
list the integrated isoscalar octupole strengths of the broad neutron mode and the 
surface vibrational mode in ^'^~^^Ni. The obtained isoscalar octupole strength for the 
broad neutron mode is 5.85 x 10^ fm^ in ^^Ni. Note that the energy-weighted sum 
of this strength is 1.7 percent of the classical isoscalar octupole sum-rule value^^^'^^-* 
^cl^ ~ ("^ (^^)n ~'~ ^ ■ '^^^ isoscalar strength of the surface vibrational 

mode is 6.74 x 10^ fm^ in the same nucleus, and it exhausts 29 percent of the classical 
sum rule value. The strength of the broad neutron mode is smaller than that of the 
surface vibrational mode by about a factor of ten. 

3.2. Transition densities of low-lying octupole modes in ^"^M 

The natures of the sharp-peak mode and the broad neutron mode become more 
evident by looking into transition densities. We here analyze three kinds of transition 
densities 

P%ir) = {^,\Y,i^l{ra)i^,{ra) |^o) = l^Af (3-1) 



P^^(r) = {<P,\ ^l:l{r T)V^|(r j) |^o) = Yl,,{r)P^^l{r), (3-2) 
P,fir) = {<P,\ i^,{r i)i^,{T T) l^o) = YlM{r)PtH{r), (3-3) 



where plq{r) is the usual particle-hole transition density while Pfq{r) and P^^{r) 
are the transition densities for particle-pair and hole-pair, respectively, for either 
neutrons or protons {q = n,p). They are evaluated as 

C 

P%L{r) = ^lm6pqL{r,uJi), (3-4) 

C 

Pinii^) = ir—2^^i^P+,gLir,uJi) - 6p-^qLir,uji)), (3-5) 



Pltiir) = -^lui{Sp+,qL{r,uJi) + 6p.,gL{r,iOi)) (3-6) 



,hh f^-^ ^ _^ 
27rr" 

in terms of the solutions of the linear response equation ()2-2p . Here C is a normaliza- 
tion constant, which is fixed so that the transition amplitude Migi = J (r) pj^^{r)r'^ dr 
gives the integrated isoscalar octupole strength of a mode under consideration by 
the standard definition S(IS3) = 7| Afig^^p. 

In Fig. \5\ we show the transition densities associated with the surface vibrational 
mode. We evaluate them at E = 4.15 MeV, approximately at the peak energy. It 
is seen in Fig. [5] that the particle-hole transition densities Piqi{r) of both neutrons 
and protons exhibit large and in-phase amplitudes at around the nuclear surface 
(the calculated matter r.m.s. radius of ^''Ni is 4.28 fm). We thus confirm that the 
mode is typical of the surface vibration in which neutrons and protons give coherent 
contribution. It is seen also that the amplitude of the particle-hole transition den- 
sity pfgi^ir) is significantly larger than those of the particle-pair transition density 
P[gl{r) and the hole-pair transition density P^^^ir) both for neutrons and protons. 
{PfqL^'^) ~ PiqLi'^) ~ ^ protous is a trivial consequence of the zero proton pairing 
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gap Ap = 0.) The dominance of the particle-hole amplitude is consistent with the 
observation that the mode is generated by the particle-hole residual interaction (cf. 
SecEU. 

Figure [6] is the transition densities of the broad neutron mode, evaluated at E = 
3.0 MeV. They have characters different from those of the surface vibrational mode. 
We observe here two distinct features. Firstly, the particle-pair transition density 
Pf^l{r) of neutrons is significantly larger than the neutron particle- hole transition 
density pfgi^ir) in the nuclear exterior region. The ratio between the two transition 
densities is approximately a factor of 3 at r = 10 fm, and a factor of ~ 10 at 
r = 15 fm. This indicates that the mode is characterized, especially in the external 
region, by motion of neutron pairs rather than by particle-hole excitations. In other 
words, the neutron pair correlation is the main character of this mode. Secondly, 
both the particle-hole and particle-pair transition densities of neutrons exhibit a very 
long tail extending to r ~> 20 fm, especially in the particle-pair transition density. 
This indicates that neutrons in weakly bound and continuum orbits participate in 
forming this mode. In addition to these two features, we observe also that the proton 
amplitude of the particle-hole transition density is considerably smaller than that of 
neutrons. This is in accordance with the dominance of the neutron strength already 
observed in the previous subsection. 

We would like to emphasize that the residual pairing interaction playing role 
in the QRPA equation brings the correlation to this mode, and that the pairing 
mean- field alone is not sufficient. This is seen in the difference between the solid and 
dashed lines in Fig. \7\ where the dynamical pairing effect, i.e., the RPA correlation 
due to the residual pairing interaction, is either included or neglected. It is seen that 
the dynamical pairing effect enhances the particle-pair transition density by a factor 
of ~ 2. This indicates significant configuration mixing effect originating from the 
residual pairing interaction. The important role of the pair correlation is seen even 
when the particle- hole residual interaction is neglected (the dotted curve in Fig. [Tj) . 

The large dynamical pairing effect is analyzed in more details. In Fig. [8] we 
show how the transition densities of the neutron mode changes if we include only 
a part of neutron single-particle orbits with lower orbital angular momenta i.e., 
including only orbits up to a smaller cut-off Icut on I. This is the same analysis 
that we performed for the soft dipole mode.^^-* It is seen that the convergence of 
the particle-pair transition density Pfgx.i''^) neutrons with respect to / is slow, and 
orbits with large / contribute coherently to produce the dynamical pairing effect. 
There is sizable contribution even around / ~ Imax = 17. Note that the precise 
treatment of the continuum states, guaranteed in the continuum QRPA approach, is 
essential to describe the correlation since high-/ orbits with I > 4 are all continuum 
states. 

Based on the above two features, i.e., the large effect of the dynamical pairing 
and the high-/ contribution, we can argue that a neutron pair moving in the neutron 
mode exhibits a spatial correlation at small relative distance between the paired 
neutrons. The configuration mixing involving high-/ orbits to a certain value, say 
Icorr, means a two-particle correlation with a small opening angle Ocorr ~ ^/Icorr-, and 
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hence a large Icorr means a spatial correlation at small relative distance. These are 
the same features seen in the case of the soft dipole excitation/^^ and we conclude 
that the di-neutron correlation appears also in the neutron mode in the octupole 
response. 
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Table I. The isoscalar and electric octupole strengths i3(IS3) and B(E3) of the neutron mode and 
the surface vibrational mode. The energy weighted sum 5^^(183) of the isoscalar octupole 
strength in unit of the classical sum-rule value Sf^. is also listed. The neutron vs proton ratio 
Mn/Mp of the transition amplitudes as well as the peak energy is also shown for the surface 
vibrational mode. The values of as well as the one- and two-neutron separation energies 
Sin and S2n are also tabulated. 



3.3. Isotopic dependence of the neutron mode 

The strength of the neutron mode increases significantly as the system ap- 
proaches toward the neutron drip line. Figure [9] shows the strength functions cal- 
culated for 80'82,84,86]\jj^ The calculation is the same as that in §3.1[ The calculated 
one-neutron separation energy in these isotopes are Ein = 2.79, 2.54, 1.86, and 1.32 
MeV in the corresponding order, and shown in FigJO] with an arrow. It is seen in 
all the nuclei that there exists above Ein a broad distribution of predominant neu- 
tron strength, which corresponds to the neutron mode. Clearly the magnitude of 
the strength increases monotonically with increasing N. This is also seen in the 
integrated isoscalar strength i?(IS3) of this mode, listed in Table HI 

The transition densities of the neutron mode, evaluated at E = 4.0, 3.5, 3.0, 
and 2.5 MeV for 80,82,84,86]vgj^ respectively, are shown in Fig. [TOj In the top panels we 
see that the tail of the particle-hole transition density grows as the neutron number 
increases. It is known that such a long tail can be realized for particle-hole transitions 
from weakly bound orbits, whose wave functions have long tail, to continuum orbits 
with small kinetic energies. The long tail enhances the octupole strength as the 
octupole operator has the radial form factor oc giving a heavy weight at larger 
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distances. We thus see that the increase of strength with increasing N originates 
from the effect of weak binding of neutrons. 

Looking at the particle-pair transition densities shown in the bottom panels in 
Fig. [TOl we observe that the particle-pair transition density varies with more 
drastically than the particle-hole transition density. In '^^Ni and also in ^'^Ni this 
transition density does not show exponential decay in the outside of the nucleus 
r > 10 fm, but it rather shows an oscillatory behaviour with its maximum amplitude 
at around r ~ 15 fm far outside the nucleus. This suggests significant emission of a 
neutron-pair from the nucleus. This is of course related to the feature that the two- 
neutron separation energy E2n decreases from E2n = 3.84 MeV in ^'^Ni to a small 
value 1.37 MeV in ^^Ni. The di-neutron correlation in the neutron mode becomes 
more significant as we approach the neutron drip-line. 

3.4. Model dependence 

Let us now investigate how our predictions on the neutron mode depends on the 
model parameters. For this purpose, we shall compare the above results with those 
obtained with another Skyrme parameter set SkM* and those with a model adopt- 
ing a Woods-Saxon potential^^^ instead of the Skyrme-HFB self-consistent mean- 
fields. When calculating the octupole response for these models, we use the same 
mixed-type DDDI, but the interaction parameter Vq is adjusted separately to re- 
produce the average pairing gap of neutrons in stable nuclei as in the SLy4 case: 
Vq = -265 and -280MeV fm'^ for SkM* and WS, respectively. The obtained aver- 
age neutron pairing gap in ^'^Ni is 1.349, 0.797 and 0.569 MeV for the Woods-Saxon, 
SLy4 and SkM* models, respectively. 

The octupole strength function in ^^Ni obtained for the SkM* and Woods-Saxon 
models are shown and compared with that for SLy4 in Figllli It is seen that the 
neutron mode depends rather sensitively on the models. The Woods-Saxon model 
produces a significantly larger strength than SLy4 while it is smallest for SkM* 
among the three models. Figure [T2] is a comparison of the transition densities of the 
soft neutron mode of the Woods-Saxon model and that of the SLy4, both evaluated 
at ii^ = 3.0 MeV. (The transition densities for SkM* is not shown here as the strength 
of the neutron mode itself is weak.) The Woods-Saxon model exhibits enhanced am- 
plitudes of the transition densities, especially of the particle-pair transition density, 
compared with the SLy4 model. 

We can relate the model dependence to differences in the single-particle lev- 
els and the Fermi energy of neutrons. The neutron Fermi energy in *^Ni is = 
-0.722, -1.183 and -2.338 MeV for the Woods-Saxon, SLy4 and SkM* models, 
respectively. The neutron single-particle levels are shown in Fig|13l Clearly the 
last neutrons are less (more) bound in the case of the Woods-Saxon model (the 
SkM* model). Thus the model dependence of the low-lying octupole strength can 
be explained in terms of the weak-binding effect which increases the strength of the 
neutron mode. 

We can speculate that the neutron mode depends also on the effective pair in- 
teraction. Let us examine explicitly dependence on the effective pair interaction. 
For this purpose we compare results obtained with the mixed type DDDI with an- 
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other calculation using the density- independent delta interaction, which is defined 
by Eq. (j2-ip but with r] = 0. (The same SLy4 is used as the Skyrme parameter set.) 
The strength Vq = — 215MeV fm^^ is adjusted in the same way as for the mixed 
type DDDI, i.e. to reproduce the average neutron paring gap in stable Ni isotopes. 
The density independent delta interaction is often called the volume pairing. The 
volume pairing produces average neutron pairing gap of 0.467 MeV in ^"^Ni, slightly 
smaller than that in the mix-type DDDI. As far as the octupole strength function is 
concerned, there is no significant dependence on the type of the pairing interaction 
as we see from the comparison in Figllll There is, however, large difference in the 
particle-pair transition density Pfqi^ir) of neutrons, as shown in Fig. [TH The am- 
plitude of the particle-pair transition density Pfqi^{r) in the volume pairing case is 
smaller by a factor of ~ 2 than that of the mixed-type DDDI. It is seen also that 
the dynamical paring effect is significantly smaller in the volume pairing case. These 
results indicate that the di-neutron correlation in the neutron mode is sensitive to 
the choice of the effective paring interaction: effective pairing interactions such as 
the mixed-type DDDI that has stronger interaction strength at low densities give 
stronger di-neutron correlation in the neutron mode. This feature is in agreement 
with our previous finding for the soft dipole excitation. 

3.5. Comparison with soft dipole excitation 

The analyses in the preceding subsections revealed that the weak binding effect 
and the di-neutron correlation in the octupole neutron mode are similar to those seen 
in the soft dipole excitations. ^^-^ In this subsection let us make a more explicit and 
quantitative comparison between the dipole and octupole cases. For this purpose 
we calculate the electric dipole response using the same Skyrme-HFB + continuum 
QRPA model. (The calculations in our previous work^^-* is based on the Woods- 
Saxon model, and not suitable for direct comparison with the present calculations.) 

The calculated electric dipole strength function dB(El)/dE is shown in FigllSi 
In the same energy region as the octupole neutron mode (i.e. just above the one- 
neutron separation energy), there emerges a broad distribution of the El strength 
corresponding to the soft dipole excitation. It is seen also that the monotonic increase 
of the strength with N is in parallel with that of the octupole neutron mode. If we 
evaluate the B(E1) strength integrating within an energy interval [Ein, i?i„ + 5MeV] 
above the one-neutron separation energy (as done in Ref.^^-*), the obtained El 
strength is B{El;Ogs 1") = 0.83,1.59,2.23,3.33 e^fm^ for 80,82,84,86]s^i^ respec- 
tively. The energy weighted sum in the same energy interval is 1.9, 3.3, 4.0 4.8% 
of the classical (TRK) sum rule value for 80.82,84,86]\j-^_ The El strength and the 
energy weighted sum in a smaller energy interval [Ein, Ein + 1.5MeV] are 0.11, 0.24, 
0.35, 0.87 e^fm^ and 0.15, 0.30, 0.36 0.68% for the same isotopes. If we compare 
the strengths of the octupole neutron mode and the soft dipole mode in terms of 
the fraction of the energy weighted strength to the classical sum rule, 1.7% for the 
octupole neutron mode and 4.0% for the soft dipole in '^''Ni (0.4% if the same energy 
interval [Ein, Ein + 1.5MeV] is used) are comparable. 

In Fig. [TBI we show the transition densities of the soft dipole mode. The tran- 
sition densities are evaluated at the same excitation energy E = 3.0 MeV and in 
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the same energy interval [Ein,Ein + 1.5 MeV] as is done for the octupole neutron 
mode. Comparing the transition densities of the octupole (FiglH]) and dipole (Fig ll6p 
modes, both are similar in that the particle-pair transition density P^^iir) is most 
dominant, and also that the dynamical pairing effect and coherent contributions of 
the high-/ orbits play significant role to enhance the particle-pair transition density. 
If we compare the absolute magnitudes of the transition densities, it is seen that the 
magnitude of the particle-pair and particle-hole transition densities of the octupole 
neutron mode is smaller than those of the soft dipole mode by a factor of ~ 2 at 
around r = 10 fm where the particle-pair amplitude is largest. Apart from this 
difference, the significance of the di-neutron correlation is comparable in both cases. 

The most noticeable difference between the octupole and dipole response is that 
in the octupole response the neutron mode overlaps with the surface vibrational 
mode present in the same energy region, and the strength is overwhelmed by the 
latter. In contrast, the soft dipole excitation in the dipole response is well separated 
from the other mode of excitation, the giant dipole resonance. This may make it 
more difficult to identify experimentally the octupole neutron mode. 

3.6. Surface vibrational mode 

Let us mention briefly the octupole surface vibrational mode, corresponding to 
the sharp peaks around E ~ 5.6 — 3.8 MeV in ^''"'^^Ni. The peak energies as well 
as the isoscalar and electric octupole strengths associated with the surface vibra- 
tional mode are listed in Table HI The most noticeable feature is that the isoscalar 
strength increases steeply with increasing from 80 to 86 approximately by a factor 
of three. We also see that the neutron vs. proton ratio M„/Mp of the transition 
amplitudes increases more steeply than the the nominal ratio N/Z. The double 
ratio [Mn/Mp] /{N / Z) increases from 1.41 to 1.88 with increasing N. These indi- 
cate that the enhanced collectivity of the surface vibrational mode is due mainly 
to neutron contributions. Note that the ratio of the energy weighted sum of the 
isoscalar strength to the classical sum rule value Sfj^ stays constant around 25-33 
%. This is partly because the increase of the strength is compensated by the decrease 
of the excitation energy (see Table H]) , and partly because the sum rule value S^i'^ 
itself increases with N due to increasing radial expectation value of neutrons. 

This suggests that the increase of the isoscalar and neutron strength of the surface 
vibrational mode is regarded as a kind of softening caused by the weak binding of 
neutrons. 

The enhanced collectivity of the octupole vibrational mode in neutron-rich nuclei 
close to the drip-line is pointed out in Ref.,^^-* which however analyzed a doubly- 
closed-shell nucleus ^''Ca using the Skyrme-HF plus continuum RPA without the 
pair correlation. Our results suggest that the enhanced collectivity is generally seen 
in nuclei near the neutron drip line. A large deviation of the Mn/Mp ratio from the 
nominal ratio N/Z is pointed out in a Skyrme-BCS+ QRPA calculation^^) for the 
neutron-rich oxygen isotope ^^O, but with much smaller deviation in less neutron- 
rich isotopes ^^~^^0. We refer also to Ref.,^''-' which discuss enhanced collectivity 
due to the weak binding effects in the case of the low-lying quadrupole vibrational 
mode. 
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§4. Conclusions 

We have investigated the low-lying octupole exeitations of the neutron-rieh Ni 
isotopes beyond the = 50 shell closure using the eontinuum QRPA based on the 
Skyrme Hartree-Fock-Bogoliubov mean-fields. In addition to the surface vibrational 
mode of the l-Hco character, a broad strength distribution of predominantly neutron 
component appears just above the neutron separation energy. This broad neutron 
mode exhibits the following distinctive features, (i) The transition amplitude for 
neutron pair density is significantly larger than that of the usual particle-hole transi- 
tion density of neutrons. The neutron pair correlation is therefore the most essential 
aspect characterizing the neutron mode, (ii) The large neutron-pair transition den- 
sity originates from a coherent contribution of the neutron high-/ orbits, which points 
to the presence of the spatial correlation among neutrons involved in this mode, (iii) 
The transition densities of the neutron mode display a long tail extending in the 
outside of the nucleus. The strength of this mode increases monotonically with in- 
creasing from 52 (^'^Ni) to 54, 56 • • • . Both indicate that the mode originates 
from the weak binding of neutrons. Prom these features we conclude that the spatial 
di-neutron correlation shows up also in the neutron mode in nuclei near the neu- 
tron drip-line, similarly to the soft dipole mode. This supports our expectation that 
soft modes having the di-neutron character may emerge generically in medium mass 
neutron-rich nuclei near the drip-line. It is interesting to investigate this generality 
in more details, e.g. by looking into modes with other multipolarities. Such an 
analysis is in progress, and will be reported elsewhere. 
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Mass Number 

;. 1. The average neutron pairing gap of nickel isotopes. The solid line shows the values obtained 
in the present HFB calculation using SLy4 and the mixed-type DDDI. The dashed line is the 
experimental value extracted from the nuclear masses using the three-point formula. 
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Fig. 2. The isoscalar octupole strength dB{lS3)/dE, plotted with the solid curve, obtained for 
**Ni using SLy4 and the mixed-type DDDI. The smoothing parameter is e = 0.2 MeV. The 
unperturbed strength is also shown with the dotted curve. The arrows indicate the one- and 
two-neutron separation energies. 
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Fig. 3. The low-energy part of calculated octupole strengths in Ni. The isoscalar, neutron and 
proton strengths are plotted with the solid, dashed and dotted curves, respectively. The smooth- 
ing parameter is e = 0.05 MeV. See also the caption of FigO 
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Fig. 4. The effects of the residual interactions on the low-energy part of the isoscalar octupole 
strength in **''Ni, shown in Fig. [S] The solid curve is the full calculation while the dashed 
and dotted curves are the results neglecting the residual pairing and particle-hole interactions, 
respectively. The dot-dashed curve is the unperturbed strength where both of the residual 
interactions are neglected. 
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Fig. 5. The transition densities evaluated at i? = 4.15 MeV for tlie octupole response of ^''Ni (cf. 
Fig. m. In the upper panel shown is the particle-hole transition density r"^p'qi(r) weighted 
with the volume element r^. The neutron and proton amplitudes are plotted with the solid 
and dashed curves, respectively. The particle-pair transition density r'^Pl^i^ir) and the hole- 
pair transition density r^Ptq'hi'r) of neutrons are shown in the bottom panel with the solid and 
dotted curves, respectively. P[qi^{r) and -Pigi,(r) for protons (the dashed curves) are exactly 
zero. 
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Fig. 6. The same as Fig. [5] but for the transition densities evaluated at i? = 3.0 MeV. 
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E — 3.00 MeV. Here the effects of the residual interactions are shown. The dotted curve is the 
case where the residual particle-hole interaction is neglected while the dashed curve is the case 
where the dynamical pairing effect, i.e. the residual pairing interaction, is neglected. The full 
calculation is also shown with the solid curve for comparison. See also the caption of Fig. [5] 
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Fig. 8. The dependence of the neutron transition densities of the neutron mode, evaluated at _E = 
3.0 MeV, on different values of Icut ~ 7, 9, ■ ■ ■ , 17ft. The upper and lower panels are for the 
particle-hole and particle-pair transition densities, piqL{f) and r^Piqi^{r), respectively. The 
transition densities obtained neglecting the dynamical pairing effect are also plotted with the 
thin solid curve for comparison. 
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Fig. 9. The octupole strength functions obtained for *®Ni. See also the caption of Fig. [3] 
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Fig. 10. The neutron transition densities of the neutron mode obtained for 



^Ni. The upper 



and lower panels are for the particle-hole and particle-pair transition densities, r'^p*J^(r) and 
r'^Pfqi^{r), respectively. We calculate the transition densities in each nucleus at the excitation 
energy listed in the Figure. The transition densities obtained neglecting the dynamical pairing 
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Fig. 11. The octupole strength functions in Ni, obtained with the Woods-Saxon model (left 
top) and the Skyrme HFB model with SkM* (right top). The mixed- type DDDI is used in 
both cases. The bottom panels show the strength functions for the Skyrme parameter set SLy4 
but with different pairing interaction: the mixed-type DDDI (bottom left) and the volume-type 
density-independent delta interaction (bottom right). See the caption of Fig.[3]for the definition 
of curves. 
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Fig. 12. The comparison of the neutron transition densities P%l{''') and r'^Pfqi^{r) of the neutron 
mode in **Ni obtained with the Woods-Saxon model (the thick curves) vs. those with the Skyrme 
HFB with SLy4 (the thin curves). The mixed- type DDDI is used and the excitation energy is 
E = 3.0 MeV in both cases. The dashed curves in the bottom panel represent the results 
obtained neglecting the dynamical pairing effect. 




Fig. 13. The neutron single-particle energies for Ni in the Woods-Saxon model and the Skyrme 
Hartree-Fock model with SLy4 and SkM* . The dashed line indicates the neutron Fermi energy. 
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Fig. 14. The dependence of the neutron transition densities r^p*Jj^(r) and r^Pfq^i''') of the neutron 
mode on the pairing interaction. Wo compare results obtained with the mixed-type DDDI (the 
thin curves) and the volume- type density-independent delta interaction (the thick curves) . The 
excitation energy \s E = 3.0 MeV in both cases. The Skyrme parameter set SLy4 is used. 




Fig. 15. The electric dipole strength function dB{Ei'i) / dE in **'~*^Ni, shown in the top panels, 
obtained with the Skyrme parameter set SLy4 and the mixed-type DDDI, compared with the 
isoscalar octupole strength function dB{IS3)/dE shown in the bottom panels. The smoothing 
parameter is e = 0.2 MeV. 
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16. The same as Fig|8] but for the soft dipole excitation. Here the angular momentum cut- 
off is Icut = 5, 7, 9, • ■ ■ 17ft. The transition densities are evaluated ai E = 3.0 MeV using the 
mixed-type DDDI and SLy4. 



